%
% IVP que se debe resolver (el del sistema a lazo abierto)
% Parametros de entrada:
%
%   t:      instante de tiempo actual
%   x:      vector con los valores actuales de x1 y x2
%   params: vector de parametros del sistema
%   u:      nombre de la funcion que modela al controlador
%   e:      ?
%   ts:     vector de tiempos hasta el momento
%
% Salida: un vector con x1' y x2'
%

function rta = ivp(t, x, params, u, e, ts)

b = params(1);
i = params(3);

x1 = x(1);
x2 = x(2);

ut = feval(u, t, params, e, ts);

x1_prima = x2;
x2_prima = (-b*x2+ ut)/i;

rta = [x1_prima x2_prima];

